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Abstract 

We present first realistic numerical simulations of 3D 
radiative convection in the surface layers of main sequence 
A-type stars with T eS = 8000 K and 8500 K, logg = 4.4 
and 4.0, recently performed with the C0 5 BOLD radia- 
tion hydrodynamics code. The resulting models are used 
to investigate the structure of the H+Hel and the Hell 
convection zones in comparison with the predictions of lo- 
cal and non-local convection theories, and to determine 
the amount of 'overshoot' into the stable layers below 
the Hell convection zone. The simulations also predict 
how the topology of the photospheric granulation pattern 
changes from solar to A-type star convection. The influ- 
ence of the photospheric temperature fluctuations and ve- 
locity fields on the shape of spectral lines is demonstrated 
by computing synthetic line profiles and line bisectors for 
some representative examples, allowing us to confront the 
3D model results with observations. 

Key words: Stars: A-type - Stars: convection - hydrody- 
namics - radiative transfer 



1. Introduction 

In comparison with the Sun, convection in the envelopes of 
A-type stars is a rather inefficient energy transport mech- 
anism. According to local mixing- length theory (MLT, 
Boh m-Vitense 1958|) convection is confined to two sepa- 
rate shallow convection zones near the stellar surface. Un- 
fortunately, the structure and convective efficiency of these 
layers depends sensitively on the choice of the -unknown- 
mixing- length parameter. Another problem with MLT is 
that it cannot describe convective overshoot. To overcome 
these difficulties, 'parameter-free', non-local convection the- 
ories have been developed and applied to A-type stars 
QKupka fc Montgomery 2002| KM02). 

Radiation hydrodynamics simulations of stellar con- 
vection constitute an independent approach. Up to now, 
realistic simulations of surface convection in A-type stars 
have been restricted to 2D preytag et al. 1996| ). 3D sim- 
ulations are challenging, because the short radiative time 
scales in the atmospheres of these stars enforce an exceed- 
ingly small numerical time step, and hence make convec- 
tion simulations for A-type stars much more time con- 



suming than for the Sun. On the other hand, A-type stars 
have the advantage that the entire convective part of the 
envelope can be included in a single simulation box with 
a simple closed lower boundary. 

In the following, we present first results of 3D hydrody- 
namical convection simulations for main-sequence A-type 
stars (T cff = 8000 K and 8500 K, log 5 = 4.4 and 4.0, solar 
metallicity), and compare them with the aforementioned 
convection theories. We also address the interesting ques- 
tion of whether the hydrodynamical models can reproduce 
the peculiar line profiles and lines asymmetries (inverse bi- 
sector C-shape) observed for slowly rotating A-type stars 
with T cff around 8000 K IjLandstreet 1998|l . 

2. 3D Hydrodynamical Convection Models 

The numerical simulations presented here were performed 
with C0 5 BOLD, a 3D radiation hydrodynamics code de- 
signed to model stellar convection (see |Freytag et al~2 002 
or Wcdcmey er et al. 2004| for details). The integration of 
the equations of hydrodynamics is based on a conservative 
finite volume approach using an approximate Riemann 
solver of Roe type together with a van Leer reconstruction 
scheme. The Roe solver was modified to handle an exter- 
nal gravity field and an arbitrary tabulated equation of 
state (EOS). For the present application, we use a realis- 
tic EOS table accounting for partial ionization of hydrogen 
and helium (as well as H2 molecule formation). 

The 3D non-local radiative transfer is solved on a sys- 
tem of long rays, employing a modified Feautrier scheme. 
Using a realistic Phoenix-OPAL Rosseland mean opacity 
table, we adopt the grey approximation in this exploratory 
study. Strict LTE is assumed (no scattering), and radia- 
tion pressure is ignored. 

The simulations are performed on a Cartesian grid 
with variable cell size in the vertical direction. We ap- 
ply periodic lateral boundary conditions, while top and 
bottom boundaries are 'closed'. 

3. Results 

The results shown here are derived from convection sim- 
ulations for stellar parameters T e g = 8000 K, log g = 4.4 
(model 1), and T ofr = 8500 K, log 5 = 4.4 (model 2), re- 
spectively. A representative snapshot from each of these 
sequences is displayed in Fig. ^ Similar calculations have 



Proc. 13th Cool Stars Workshop, Hamburg, 5-9 July 2004, F. Favata et al. eds. 



2 



MODEL. aeOa44n10_3Db:.enJ T ; „-S™C I- 4.40D 

IfeD □ .25 1300 Mm y^H DDDDD □.□□□□0 Mm fc11 253815 D729S6351 Mm 
ij j k = 180 j 1BD x90 EOS: fDS_mmDDJeos V: phuenix_opal_grey opta 



Specific Entropy 




MODEL: aB5g44n15_12end T^5DD K bgfg> 4.400 
X^O Q70D0DQ 25 1300 Mm y^D D7DDDDB . Q D70000B Mm 4^12 774QS3 3 35621 54 Mm 
xjxk = 130 x IBOx 110 EOS:eo5_mm00J5eos Kl: phoenK_opal_grey opta 

3DA-type star at85g44nl 5. Time=21603s: Temperature 




Figure 1. Arbitrary snapshots from two 3D simulations of 
convection in the surface layers of A-type stars, showing 
the velocity field and temperature distribution in a vertical 
slice. Top: MODEL 1: T cS = 8000 K, \ogg = 4.4, geomet- 
rical size 25.2 x 25.2 x 12.0 Mm, 180 x 180 x 90 grid cells, 
vertical optical depth range —3 < log tr OS s < 4, covering 
« 10 pressure scale heights. Bottom: MODEL 2: T c jf = 
8500 K, log g = 4.4, geometrical size 25.2 x 25.2 x 15.0 Mm 
, 180 x 180 x 110 grid cells, —6 < logTR OSS < 4, covering 
w 16 H v . 



been performed for log g = 4.0, but arc not shown here. 



Further details can be found in Freytag & Steffen (2004) 
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Figure 2. Mean specific entropy (s(z)) x ,y,t for Model 1 
(top) and Model 2 (bottom). In both cases, the en- 
tropy stratification indicates two separate convection zones 
(ds/dz < 0, shaded). The borders of convective instabil- 
ity according to local MLT models are indicated by thin 
(a=0.5) and thick (a =1.0) arrows. 



tions, although considerable differences remain, especially 
in the energy fluxes. In contrast to the findings by KM02, 
our 3D models do not show anywhere a positive kinetic en- 
ergy flux. Tabled lists a some key numbers characterizing 
the different kinds of models. 

Fig. 0] demonstrates that 'overshoot' below the Hell 
convection zone is substantial. The exponential tail of the 
velocity field is clearly seen in model 2, where the veloc- 
ity scale height in terms of the pressure scale height at 



3.1. Vertical Structure, comparison with MLT the bottom of the Hell CZ is H v /H p w 0.4. Model 1 is 

not deep enough to include the exponential part of the 
overshoot region. We estimate H v /H p <, 0.7. 



In Figs. [3- El the mean vertical structure of the 3D hydro- 
dynamical simulations, obtained by horizontal and tem- 
poral averaging, is compared with the results of standard 
MLT, in the version described by |Mihalas (1978)1 for 

mix- 
ing length parameters a=0.5 and 1.0. Clearly, MLT does 
not even approximately match the hydrodynamical results, 
no matter what value of a is chosen. We note that the non- 



local convection models by Kupka & Montgomery (2002) 
are qualitatively more similar to the hydrodynamical solu- 



3.2. Horizontal Structure 

The horizontal structure emerging from our 3D simula- 
tions of surface convection in A-type stars is evident from 
the intensity images displayed in Fig. 03 Obviously, the 
flow topology is qualitatively similar to that of the so- 
lar granulation, i.e. isolated hot up-flows (granules) are 
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Figure 3. Mean vertical velocity y/ (V^(z)) x ,y,t) compared 
with MLT models. According to 3D hydrodynamics, both 
convection zones are connected by overshooting flows. The 
extended exponentially decaying flow field below the He II 
convection zone is also a result of 'overshoot'. 



Figure^. Mean enthalpy flux (pV z h} x ^ y j, compared 
with mixing-length results. Like in the solar granula- 
tion, the 3D simulations give a mean kinetic energy flux 
(pV z V 2 /2) XtVt t which is directed downwards at all heights. 



separated by a network of connected cool down-flows (in- 
tergranular lanes). Due to a more efficient radiative en- 
ergy exchange, the granules at the surface of A-type stars 
are relatively larger than on the Sun, and seem to show 
less sub-structure. In fact, the filling factor of dark areas 
(where / < I) is fd « 0.34 for both model 1 and model 2, 
compared to fd ~ 0.53 for the Sun. 

3.3. Synthetic Line Profiles 

For the snapshots shown in Fig. [5] we have computed syn- 
thetic line profiles both for vertical rays (disk-center) and 
for integrated light (flux) under the assumption of LTE. 
The resulting disk-center line profiles and line bisectors of 
Fel A 6265.13 A are presented in Fig. |SJ The considerable 
photospheric velocities and temperature fluctuations in- 
duce a distinct asymmetry of the emergent line profiles. 
We have investigated several snapshots and different spec- 
tral lines, and found that the line bisector always exhibits 
a solar-like C-shape, but with a larger excursion to the 
red near the continuum. In the flux spectra, the line bi- 



sectors span typically 2 km/s and 1 km/s for models 1 and 
2, respectively. This is of the same order of magnitude as 
observed by Landstreet (1998)) but the asymmetry is in 
the opposite direction. The C-shape persists for the sim- 
ulations with log g= 4.0. 



4. Conclusions 

The analysis of our 3D hydrodynamical simulations indi- 
cates a severe failure of the standard local mixing-length 
theory in the regime of A-type star shallow surface convec- 
tion. The non-local convection model by Kupka & Montgomery (2( 
gives a much better description of the velocity field, but 
alarming differences remain in the energy fluxes. Over- 
shoot below the He II convection zone is found to be sub- 
stantial, in basic agreement with KM02. 

According to the simulations, the granulation pattern 
forming at the surface of A-type stars has a solar-like flow 
topology, with granules that are relatively larger than on 
the Sun and seem to show less sub-structure. Synthetic 
LTE line profiles based on the current 3D convection mod- 
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Continuum surface intensity 0- 5265 A) 



Con'tlnulim au.rfaGe intensity [X 6265 A) 





Continuum surface intensity (>. 6265 A) 




Figured. Emergent eontinuum intensity at A 6265^4 resulting from 3D hydrodynamical simulations of the solar gran- 
ulation (left, 5/ rms =14.6%) and of surface convection in main-sequence A-type stars with T q r = 8000 K (middle, 
« rms =12.7%) and T cff - 8500 K (right, 5J rms =7.9%). 
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Figure 6. Spatially resolved and averaged line profiles (top) and line bisectors (liottomj of Fe IX 6265.13 A, computed 
from the snapshots shown in Fig. 5 for vertical lines- of- sight (fj,— l). For the Sun, the gf -value has been reduced by a 
factor 100. The asymmetry of the flux profiles (not shown) is qualitatively similar. 



els of A-type stellar atmospheres show a depressed red 
wing, in apparent contradiction to the observed line asym- 
metry (|Landstreet 1998J1 . A possible reason for this dis- 
crepancy might be missing physics in the simulations (e.g. 
magnetic fields). On the other hand, we note that the 
slowly rotating A-type stars with T c g ~ 8000 K observed 
by Landstreet (1998) are classified as Am type and known 
to be spectroscopic binaries; hence they may be peculiar. 
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Table 1. Comparison of maximum convective velocity, 
Vc max ([km/s]), maximum fraction of convective en- 
ergy flux, Fc maX7 and of kinetic energy flux (Fk max = 
max.{\Fk\) / F ) for different kinds of A-type star convec- 
tion models, (u) and (I) refer to the upper (H+Hel) and 
lower (Hell) convection zone, respectively. 





Hydrodyn. 


MLT 


MLT 


KM02 




Simulation 


a = 0.5 


a = 1.0 


Theory 


T eff = 8000 K, 


log g = 4.40: 








VCmix (u) 


5.59 


3.71 


6.17 


5.48 


Vc max (1) 




0.16 


1.05 


2.48 


fc max (u) 


0.255 


0.458 


0.806 


0.100 


Fc max (l) 


0.031 


0.000 


0.012 


0.030 


-Ffc max (u) 


0.0215 






0.0011 


Ffc max (l) 








0.0012 


T eff = 8500 K, 


log g = 4.40: 








Vc max (u) 


3.90 


1.90 


5.94 


5.29 


Vc max (1) 


1.04 


0.11 


0.82 


2.70 


fc max (u) 


0.060 


0.030 


0.430 


0.019 


Fc max (l) 


0.003 


0.000 


0.004 


0.023 


Ffcmax(u) 


0.0021 






0.0004 


Ffc max (l) 


0.0002 






0.0010 



